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<N Abstract 

In synthetic biology, gene regulatory circuits are often constructed by combining smaller circuit 

qq components. Connections between components are achieved by transcription factors acting on 

(N 

promoters. If the individual components behave as true modules and certain module interface 
conditions are satisfied, the function of the composite circuits can in principle be predicted. In 
this paper, we investigate one of the interface conditions: fan-out. We quantify the fan-out, a 

O 

concept widely used in electric engineering, to indicate the maximum number of the downstream 
inputs that an upstream output transcription factor can regulate. We show that the fan-out is 
closely related to retroactivity studied by Del Vecchio, et al. We propose an efficient operational 
method for measuring the fan-out that can be applied to various types of module interfaces. We 



> 

also show that the fan-out can be enhanced by self- inhibitory regulation on the output. We discuss 



in 

ON 

o 
o 



the potential role of the inhibitory regulations found in gene regulatory networks and protein signal 
pathways. The proposed estimation method for fanout not only provides an experimentally efficient 
way for quantifying the level of modularity in gene regulatory circuits but also helps characterize 
and design module interfaces, enabling the modular construction of gene circuits. 
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Background 



Engineering relies on modular composition that is the ability to combine functional units 
with the knowledge that the intrinsic properties of each module is unaffected to a large degree 
by the composition. In biology we are less clear on the notion of a modular component, 
or at least biology has multiple definitions depending on context. Here we will define a 
module as a self-contained functional unit whose intrinsic properties are independent of the 
surrounding milieu. This definition is similar to that used in engineering. For example, the 
intrinsic properties of a TTL NAND gate pQ is unaffected (within certain design constraints) 
when connected to other TTL logic gates. That is, a NAND gate remains a NAND gate 
no matter what it is connected to. This property allows engineers to design, predict and 
fabricate complex circuits at very low cost. The question whether such self-contained and 
functionally independent modules exist at the biological cellular network level is still an 
ongoing research problem [2]. In this paper we will be concerned with the design of modular 
synthetic components [SHE] arid avoid the question of modularity in natural complex systems. 

In the most abstract sense we can define a module as follows. Given a functional unit 
M with input / and output O, we can define a relation between the input and output as 
O = M(I). Given two functional units, M\ and M 2 , where the output of Mi serves as the 
input to M 2 , then Mi and M2 are defined as modules if the relation, 2 = M2(Mi(Ji)) is 
true. This simply means that in connecting Mi and M 2 together, M 2 has no effect on the 
functional characteristics of Mi and vice versa. 

Predictable composition is of particular interest to the synthetic biology community where 
gene circuits are "wired" together via transcription factors (TFs) and corresponding promot- 
ers (Figure 1). This mode of wiring makes the physical construction of relatively complex 
networks possible [9j [10]. However the general question of whether making a connection 
between two genetic units results in a functional whole remains. In particular a number of 
issues present themselves that include independence from the surrounding milieu (also called 
orthogonality [1]): domain matching and impedance bridging. The former describes the sit- 
uation where the operating concentration range of an output transcription factor matches 
the range of the input target. Impedance bridging, which is the main topic of this paper, is 
concerned with how much a downstream target circuit can affect the functional properties 
of an upstream unit. It is not to be confused with impedance matching which is related 
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to the maximum power transfer between two circuits. There has recently been interest in 
defining impedance bridging in genetic and protein circuits [HI H2] and a related quantity 
called retroactivity was introduced by Saez- Rodriguez et al. [13] and Del Vecchio et al. [H] 
to describe the effect of one module on another. 

In electric engineering there exist guidelines and published constraints on how many 
electrical modules can be driven from a source. For example, one rule of thumb for analog 
circuits suggests that the impedance at the input should be ten times the impedance at the 
driving circuit. In digital circuits, such as TTL circuits, manufacturers will quote the fan- 
out and fan-in for a given electrical module. The fan-out indicates how many downstream 
logic gates can be connected to a given output. Exceeding these limits will potentially cause 
signal distortion in analog circuits and circuit failure in digital circuits. 

We envision the development of similar criteria for connecting two biological modules 
together in synthetic biology, and here introduce the notion of fan-out for a genetic circuit. 
We define the fan-out of a genetic circuit as the maximum number of downstream promoters 
that can be driven from an upstream circuit signal without significant time-delay or signal 
attenuation. We will show that the fan-out can be estimated by using the autocorrelation 
[T5H20] of gene expression noise [2TH21] . During the estimation procedure the system's 
retroactivity can also be measured. This fan-out /retroactivity estimation can be applied 
under quite general module interface conditions. Although our analysis is focused on genetic 
networks, the principles apply equally to signal transduction networks. 

Results and Discussion 

Module interface process 

When two synthetic gene circuits are connected, transcription factors are used to connect 
them. The reaction processes involving the transcription factors such as transcription, trans- 
lation, degradation, and downstream-module promoter regulation, will be called module in- 
terface processes (MIPs). For example, consider the repressilator [25]. Let us choose the 
TetR repressor, one of the genes comprising the oscillator, as an output of the oscillator 
module (Fig. 1). When a downstream module has tetR-operons, the MIP includes tetR- 
transcription, translation, and TetR binding/unbinding to its specific operons located in the 
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downstream module. 



Retroactivity and mapping between a module interface process and an RC-circuit 

We will investigate a MIP by mapping it to a simple electric circuit composed of a resistor 
and a capacitor connected in series (RC circuit). This mapping becomes significantly helpful 
for understanding retroactivity [H] and quantifying fan-out. 



Isolated case 



When an upstream output does not regulate any downstream promoter, we can model the 

corresponding MIP as a simple TF translation-degradation process (see Fig. 2A and B). The 

concentration of the TF, denoted by X, changes in time by following the equation 

dX , , , . . 

a- = «(«)- 7* (1) 

with ait) the translation rate and 7 the degradation rate constant. We will show how this 
process can be related to an RC circuit, where a resistor and capacitor are connected in series 
and driven by an input voltage source Vi n (Fig. 2C). The total voltage drop across both the 
resistor and capacitor is equal to the driven voltage: Vi n = RI + V out , where I denotes the 
current flowing through the resistor, and V out the voltage drop across the capacitor. The 
current is equal to the rate of charge accumulation (Q) in the capacitor: / = dQ/dt, where 
the small increment dQ causes the change in V out in proportion to dQ: dQ = CdV out , with 
C a proportionality constant called capacitance. Thus, the current / can be expressed as 
CdVout/ dt. By substituting this into Vi n = RI + V out and dividing the resultant equation by 
RC, we obtain 

dVout V n V ou t . ^ 



dt RC RC ' 

where RC is known as the response time To of the RC-circuit [26]. By comparing Eqs. ((!]) 
and (]2]), we obtain the following correspondence: X = V out , a = Vi n /RC, and 7 = 1/RC, 
and the response time is expressed as 

To = RC = -. (3) 

7 

Thus, the TF-translation-degradation process (Fig. 2B) can be directly mapped to the RC- 
circuit (Fig. 2C). 



Connected case 



Now we consider the case where two modules are connected (see Fig. 3A). We will review 
retroactivity studied by Del Vecchio et al. [2]. The retroactivity was defined as the slow- 
down in interfacial dynamics in response to interactions with downstream components. The 
retroactivity described a number between zero and one with one being the least desirable, 
i.e. the interfacial dynamics are affected most. In their analysis, they assumed that the 
binding-unbinding process of the TF is fast enough that the process can be approximated 
to be in the quasi-steady state {k on X + k Q ff ^> 7; cf. J27H2H]). They also assumed that the 
lifetime of the bound TF is much longer than that of the unbound TFs. 

Specifically, they showed that the free TF concentration X changes in time by the fol- 
lowing equation [T4"] 

— = {l-U{X)){a- 1 X), (4) 



where 7Z{X) is the retroactivity, given by 
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, with K d the dissociation 



constant for the TF with respect to the promoter, and Pt the total number of the promoters. 
They showed that 1Z is always less than 1 and non-negative. The extra factor 1 — 1Z appears 
when compared with the isolated case, resulting in the slow-down of the dynamics. More 
precisely, the slow-down is due to the decrease in the factor placed in front of X in Eq. (|4]): 
7(1 — 11), which is related to the apparent response time: 

1 

We will consider the MIP shown in Figs. 3A and B under the same assumptions as given 
by Del Vecchio et al [Ti] . To understand the retroactivity by using an RC-circuit analogy, 
consider a circuit as shown in Fig. 3C. The total capacitance becomes the sum of the two 
capacitances: CV = C + C '. Thus, the response time becomes RCt- t = RCt- The change 
in the output voltage is governed by the same equation as in the isolated case except the 
capacitance C is replaced to Ct- 



dVout Vin Vout 



dt RCt RCt 



1- r 



V,„ V, 



out 



RC RC 



(6) 



C + C 

By comparing Eqs. Q and ([6]), we find that the retroactivity is given by the relative ratio 
of the new capacitance: 
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and that the response time r corresponds to r a (Eq. (O): 



We have shown that connecting downstream promoters in the MIP is equivalent to con- 
necting extra capacitors in parallel with an existing one in the RC-circuit. Due to these 
extra capacitors, the circuit takes a longer time to fully charge all the capacitors, resulting 
in the slow-down in the circuit response time. Biologically, the bound promoters act as a 
reservoir of potentially free TFs: Whenever there is a change in the number of the free TFs, 
the reservoir quickly buffers the change in the number of free TFs |30j . Such buffering causes 
transient dynamics at the interface to slow down. 



Response time vs. number of promoters 

We now investigate how the response time r is related to the total number of promoters Pt- 
The response time is shown to increase with P^ (see the Methods section) as 

T Ptr = R(C + P T C 1 ), (9) 

where C\ is a proportionality constant satisfying CjC\ — (1 + X/Kj) 2 . The above equa- 
tion Q can be viewed as each individual promoter contributing an extra capacitance C\ to 
the total capacitance (Fig. 3C): CV = C + PtC\. The capacitance G\ of each extra capaci- 
tor is related to a unit load onto the upstream output dynamics from a single downstream 
promoter. This is an interesting result and becomes useful for proposing an experimental 
method for estimating the fan-out. 

The linear relationship between the extra capacitance and Pt (see Fig. 3C box) does not 
come from any linearization approximation, but from the fact that each downstream pro- 
moter affects the upstream as an independent effector (reservoir or sequestrator), although 
the sequestration itself is represented by a nonlinear reaction. 



Gene circuit fan-out 



We will define a gene circuit fan-out by the maximum number of promoters in a downstream 
module that the output (transcription factor) of an upstream module can regulate without 
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altering the output dynamics significantly. To exemplify how much the upstream module 
can be affected, we consider a repressilator [23] as a module and its Tet repressors as a mod- 
ule output (Fig. 1A). When the output regulates tetR promoters located in a downstream 
module, the oscillation amplitude of the tetR expression level can be significantly changed, 
e.g., 40% decrease when the number of the promoter (Pt) is changed from to 100 (Fig. IB 
and C). Our interest is here to quantify the maximum number of the promoters (fan-out) 
that the upstream module can tolerate. 

Let us quantify the fan-out by considering again the simple MIP shown in Fig. 3A. 
We consider a frequency response between the input and output voltage, V in and V out , 
respectively (Fig. 4A). In the RC circuit, the capacitor acts as a low pass filter: The gain of 
the signal (the ratio of the oscillation amplitude of V ou t to that of Vi n ) is at the maximum 
level for low frequencies and drops significantly when the circuit no longer responds as fast 
as the input signal changes (Fig. 4B). The frequency when this happens is called cut-off 
frequency (u c ) (Fig. 4B) and corresponds to the inverse of the response time: 1/RCt [HTj . 

If the upstream module functions as a synthetic oscillator, there will be a practical upper 
limit (co max ) in the oscillator's frequency (e.g., for the repressilator, u max can be the inverse 
of the repressor lifetime = log(2)/10 min -1 ~ 4 hour -1 [25J). If u max is smaller than the 
cut-off frequency u c , the oscillator output will operate in a predictable manner and the 
output signal will be passed downstream without any significant signal loss. As the number 
of the downstream promoters increases, the total capacitance increases as shown in Eq. Q 
and the cut-off frequency (u c = 1/RCt) decreases. For the cut-off frequency to be larger 
than the maximum operational frequency ou max , the total number of the promoter must be 
smaller than a certain value, which will be called the fan-out. The fan-out denoted by F w 
is obtained where u c equals u max , i.e., u c = [R(C + PtCi)] -1 = u r , 



J max- 



c 

C~i 



(10) 



In the fan-out equation (10), there are two unknown parameters: C/Ci, and ro. These can 
be experimentally estimated by performing two independent experiments with and without 
any downstream module. In each experiment we measure the corresponding response time: 
r and Tp T (by using gene expression noise as will be presented later in the Results section). 
Thus, one of the unknowns tq can be estimated. How can we estimate the other unknown 
C jC\ from rp T ? If we know a priori the copy number of the promoters Pt, we can obtain 



the value of C/C\ from Eq. (|9]). If the promoters are placed on plasmids, the copy number 
of the plasmids can be estimated depending on what type of origin of replication is used, 
and thus the copy number of the promoters Pt can be known. By calculating t /RCi, we 
can estimate the other unknown, CjC\. 



Gene circuit fan-out in more general interfaces 

Up to now we have considered a simple MIP without feedback and where the degradation 
rate is assumed to be first-order. Here we will consider the more general case and show that 



we can use the same or similar fan-out function as Eq. (10). 



Oligomer under directed degradation and self-regulation 

Consider that a TF, composed of n monomers, is tagged for degradation and that its tran- 
scription is self-regulated as shown in Fig. 5A. We can obtain the fan-out expressed by the 



same equation ( 10 ), where Tq is the time constant in the isolated case, given by the difference 
between the elasticities (32]: l/r = e 2 — £i, where e\ = dvi(x,a)/dx and e 2 = dv 2 (x)/dx 
(refer to the Methods section). This means that the fan-out can be estimated exactly in the 
same way as in the monomer case as shown in Fig. 3A. All the above results apply for the 
case that the intermediate reaction steps of the oligomerization and directed degradation 
are taken into account (refer to the SI). 



Multiple promoters having different affinities 

Consider the case that two different types of TF-specific promoter plasmids, having different 
affinities for the TF and different strength of the origin of replication. It is shown that the 
MIP can be mapped to an RC-circuit having two different capacitances connected in parallel 
to C as shown in Fig. 5B (see the Methods section). The fan-out of each promoter is shown 
to satisfy the following functional relationship between F\ and F 2 (refer to the Methods 
section): 

1/wmox = r„(l + F x -± + F 2 %, (11) 
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where Fi is the fan-out for promoter plasmids of the i-th kind, and Cj denotes the corre- 
sponding capacitance per plasmid. If there are N different kinds of promoter plasmids, we 
need to sum up all the capacitances in the above equation. We note that the fan-out is not 
a single number but is given by a functional relationship between F^s: We need to balance 
the number of plasmids of different kinds depending on its unit load on the retroactivity, 
i.e., Q/C. 

To obtain the fan-out function, we need to find three unknown parameters: r , Ci/C, and 
C2/C. tq can be estimated in the isolated case. Ci/C can be estimated in the case that only 
the i-th kind of promoter plasmids exists (under the assumption that the strength of each 
origin of replication is already known). These three independent experiments will suffice for 



estimating all the unknown parameters and proposing the fan-out function Eq. (11). 



Multiple operators 

Consider the case that the promoter region includes multiple operators specific to an output 
TF (e.g., Oi, O2, and O3) having different affinities (Fig.5C). Regardless the number of the 



operators, the same fan-out function as Eq. ( 10 ) is obtained (refer to the SI) 



Two output signals 

When two output TFs (X and Z) regulate a downstream promoter independently, i.e., if 
there is no overlap between the operator regions and somehow X does not interfere with the 
operator region of Z and vice versa, the fan-out corresponding to each output TF can be 
obtained. 



We have shown that the fan-out function like Eqs. (10) and (11) can be obtained for each 
of individual cases given above. Furthermore, for all the combinations of these individual 
cases the fan-out functions apply. 

Design scheme for fan-out enhancement 



How can we increase the fan-out? Based on the fan-out equations (10), there are two ways: 



increasing CjC\ or l/r . The way to increase the latter is to apply a negative feedback on the 
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translation of X (making e\ negative for the case shown in Fig. 5A, where 1/tq = t2 — ei) 
and a positive feed-forward on the degradation rate (increasing €2). Since the enhanced 
degradation and negative feedback decrease the concentration level of X, to prevent this, it 
is desirable to amplify the translation rate (which makes e\ more negative). This is exactly 
the same mechanism proposed by Del Vecchio et al. [14J to reduce retroactivity; when the 
retroactivity is small, the upstream output dynamics does not slow down significantly by 
connecting the output to the downstream module, meaning that the load from downstream 
to the upstream is small enough that many replicates of the load can be applied to the 
upstream without slowing down the output dynamics significantly. 

One of the mechanisms, inhibitory auto-regulation, is frequently found in Escherichia coli 
transcription factors regulating a set of operons, e.g., for amino-acid biosynthesis where a 
single TF may control multiple targets, likewise for flagella formation [33]. Such motifs are 
called sing-input-module motifs [33] . 

The concept of fan-out is not limited to gene regulatory circuits. In principle, as long as 
the same class of interface processes are found regardless of the type of biological systems, 
the fan-out and retroactivity concepts can be applied p~H |3l] . For example, in the eukary- 
otic MAPK pathway, doubly phosphorylated MAPK can activate a number of downstream 
proteins and transcription factors in the nucleus. This MAPK regulation can be described 
by the module interface process similar to the one shown in Fig. 5B (in this case, many 
promoter plasmids instead of the two). In the MAPK pathway, there is a negative feedback 
from MAPK to the phosphorylation of MAPKKK [TTl IT2] . The negative feedback increases 
the fan-out of the MAPK module thereby permitting MAPK to effectively regulate multiple 
targets and multiple homologous binding sites. 



How to measure the time constant r 



It is known that transcription factors show significant stochastic fluctuations [T7H2QI 1351 - 
[38] (for review articles, J2TH21]). Their correlation times have been measured by obtaining 
autocorrelations by in vivo time-lapse microscopy [TTH20"] . Recent numerical studies show 
that the correlation time is approximately equal to the response time of the deterministic 
case [39J and that it changes as a result of connecting two genetic systems [39, 40j. Therefore, 



by using the change in the correlation time, we can estimate the fan-out by using Eq. (10) 
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as well as the retroactivity by using 1Z = Tp ^ T ° (obtained from Eq. ^ by using Eqs. (j3]) 



Tp, 

and (|8 



Example: Fan-out/retroactivity estimation 

Let us consider the simple MIP shown in Fig. 3A as a model for TFs in E. coli. We consider 
that the average copy number of plasmids containing the specific promoters ranges from 
1 to 100. If we set the volume of E. coli roughly equal to 1/mi 3 , a copy number of one 
corresponds to 1 nM. As a result we will henceforth interchange the unit of nM with that 
of copy number. A simulation using the standard Gillespie method jH] was performed and 
the observed autocorrelation was fitted to an exponential function: G(At) = Aexp(—At/r) 
with r a correlation time (we conducted a linear fit in the log-scale in the ?/-axis and the 
normal scale in the x-axis) and 1/r obtained from the fitted slope. 

For experimentally reasonable parameter values, i.e., a = 20 nM hour -1 , 7 = 2 hour -1 , 
k on = 10 nM -1 hour -1 , and k Q ff = 10 hour -1 , we performed the stochastic simulations with 
and without any downstream-module promoter (Pt = 100 and 0). The concentration levels 
of the total TF was recorded 100 times over 2 hours, the autocorrelation of this signal fitted 
to an exponential function, and the response time measured • We obtained the error bar 
of the time constant from 10 independent replicates of the autocorrelation. 

When the translation rate was set to 20 nM hour -1 , we obtained tq = 0.52 ± 0.06 hour 
and rioo = 0.9 ±0.1 hour. We obtained C/Ci = 140 ± 20, by using 

C „ RC „ r 

-IT 



d RC T - RC t Pt - r 



=100 



where we used Ct — C = PtC\. From Eq. (10), we obtained the fan-out function for this 
MIP: 

F_ = 140[ ± 20] (iWi - l' 

If the upstream module is a synthetic oscillator with a maximum operating frequency u max = 
1 hour -1 , the fan-out becomes F = 130±20. This means that we can use promoter plasmids 
with low, medium, and high copy numbers without affecting the TF dynamics, if a single 
TF-specific operator site resides on a plasmid. The retroactivity can also be estimated from 
the measured values of r and Tioq: H = 0.4 ± 0.1. 
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If we reduce the translation rate by half (now, a = 10 nM hour -1 ) the free TF concen- 
tration decreases by half. As the concentration decreases, the retroactivity increases [TU ES] 
and the fan-out decreases. We estimated r = 0.52 ± 0.07 hour and r 100 = 1.75 ± 0.04 hour. 
For the same u max = 1 hour -1 we obtained the fan-out: F = 40 ± 1. This would mean that 
we could safely use only low copy number plasmids. The retroactivity is estimated to be 
0.70 ±0.05. 

Conclusions 

In this paper we have introduced the concept and quantitative measure of fan-out for 
genetic circuits which is a measure of the maximum number of promoter sites that the output 
TFs of the upstream module can regulate without significant slow-down in the kinetics of 
the output. In addition we have also proposed an efficient method to estimate the fan- 
out experimentally. We have shown that the fan-out can be enhanced by self-inhibitory 
regulation on the output. In the estimation process of the fanout, the retroactivity can also 
be estimated. This study provides a way for quantifying the level of modularity in gene 
regulatory circuits and helps characterize and design module interfaces and therefore the 
modular construction of gene circuits. 

Methods 

We obtain the time-constant Tp T and derive the fan-out functions for various cases. 
Monomer TF 

Consider the simple MIP shown in Fig. 3. Under the quasi-equilibrium in the binding- 
unbinding process, we obtain the concentration of the bound TF as = P^Xj{X + Kj). 
We will simplify this equation to, by introducing f(X) = X/(X + K d ), 

P b = P T f(X). (12) 

The time evolution of the total transcription factor (Y = X+P&) is governed by the following 
equation [55] : 

dY 

—— = a — 7 A . 
dt 

12 



We will obtain the response time constant of Y rather than X, because Y is a pure slow 
variable showing the dynamics of our interest [T^ fIE\ |3"§] . We denote the response time 
constant of Y by t Pt , and is obtained by taking the derivative on the right hand side of the 
above equation with respect to Y: 



d(a — jX) 
dY 



T 



dY 

dX : 



(13) 



with r = 7 1 . By using Y = X + P b , the above equation becomes t Pt = r [l + ^] . By 



using Eq. (12), we obtain the time constant: t Pt = r (l + f'(X)P T ). By comparing this 



with Eq. © we obtain C/d = K d (l + X/K d f. 



Oligomer TF under directed degradation and self-regulation 



We consider that the transcription factors are composed of n monomers described in Fig. 5A. 
We assume that the binding-unbinding process is in equilibrium, and we obtain Pt as 

P T X n 

The time evolution of Y(= X + nP&) is governed by 

dY 



dt 



v 1 (X)-v 2 (X). 



The response time constant of Y is obtained as 

-dlv^-v^X)] 1 - 1 



dY 



dY 
dX' 



where r = l/(e2 — ei) with t\ = dv\/dX and e 2 = dv2(X) / dX . By using Y = X + nPb, we 
obtain 



T~P T = T 



l + P T n 



df{X) 
dX 



where is defined to be C\jG by comparing the above equation with Eq. Q. 



Multiple promoters having different affinities 



We consider the case that two different types of TF-specific promoter plasmids, having 
different affinities for the TF and different strength of the origin of replication. The TF is 
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assumed to be a monomer. The concentration of the TF bound on the promoter of each type 
(i) is given as Pu = x+K d = PTifi(X) for i — 1,2. The response time constant t Pt1jPt2 is 
given by Eq. fll3~] ). By using K = X + P w +P b2 , we obtain t PtuPt2 = r [l + ^ + ^] , and 
rewrite this as t PtuPt2 = r [1 + /J(X)Pri + /2(X)P T2 ] by using P w = P Ti fi(X). Finally by 



equating t PtuPt2 to l/u max , we obtain Eq. (11). If there were n-different types of promoter 
plasmids, the above equation is changed by replacing the last two terms to the sum over all 
the n-types. 

In the RC-circuit representation, t Pti>Pt2 is given by RC P and the total capacitance C p 
is obtained as C T = C + C X P T1 + C 2 P T2 with Q = T f-(X)/R. This indicates that this MIP 
can be mapped to an RC-circuit having two different capacitances connected in parallel to 
C as shown in Fig. 5B. 
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Supplementary Information 



I. REPRESSILATOR BIOMODEL 

We obtained a model for the repressilator from the BioModels Database (model id: 
BIOMD0000000012) [12]. This model is composed of transcription and translation pro- 
cesses for respective lad, tetR, and cl. Inhibitory regulations among them are described by 
Hill functions. One of the repressors, TetR, is chosen as an output of the repressilator mod- 
ule. The TetR repressors are allowed to regulate a downstream module and its regulation 
is described by the binding-unbinding reactions between the TetR and its specific promoter 
located in the downstream module. The copy numbers of mRNAs of respective lad, tetR, 
and cl are denoted by X, Y, and Z, and the corresponding repressor molecules by PX, PY 
and PZ. The model process (per cell) is described as, 

K 2 

aot r +a t r 2 " 2 

K M +FZ , v kd mRNA X_ 



aotr+atr 2 2 



r A 1 



K M + PX \ k d mRNA Y 



K 2 

aotr+atr 2 +t Y 2 , , ~ 

K M+^ r v kd mRNA Z 
7- r 



kaX m A kdprotPX 

PX — > 0, 



hiY py kdprotPY^ 



kttZ kdprotPZ^ 



k on PY 2 P f 

2PY + P f ; P b , 

where the total number of the tetP promoter in the downstream module is given by 
the sum of the numbers of the free and bound promoters: Pf + Pf,. We have assumed 
that the transcription and translation processes for each repressor gene are identical. The 
parameter values used for simulations are listed in Table [T| We have modified the parameters 
to increase the retroactivity by reducing the expression levels of the repressor molecules 
(a tr : 29.97 -+ 2.97, K M : 40 -> 10, eff : 20 -+ 10). 
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aotr 
a tr 



0.03 min 1 Basal transcription rate 

3 — 0.03 = 2.97 min -1 Maximum transcription rate is set to 3 min -1 . 



K 



M 



10 



, , log(2) - -1 

kd m RNA — mm 

77"lJV/i T mRNA 

TmRNA 2 min 



hi 

eff 

Tprot 

kon 

k off 



effi 
10 



log(2) ■ -1 

&v ' mm 

mRNA 



log(2) ■ _l 

Tprot 

10 min 
0.166 min" 1 
1.66 min -1 

or 100 



Number of repressor molecules corresponding to 
the half maximal transcription 
Degradation rate constant of mRNA 
Half life time of mRNA 
Translation rate per mRNA 

Translation efficiency: Number of repressor molecules 
translated per mRNA during mRNA life time 
Degradation rate constant of repressor molecules 
Repressor molecule half life time 
Promoter-repressor binding constant [38j 
Promoter-repressor unbinding constant 
for a given dissociation constant = k ff/k on = 10 
Number of tetR promoters in a downstream module 



TABLE I: Parameters of the repressilator BioModel 



II. LINEAR RELATIONSHIP BETWEEN C AND P T 



In this section, we consider the case that the intermediate reaction steps of oligomerization 
and enzyme- mediated degradation are taken into account. 



A. Oligomer TF 

We consider a dimer TF (the results obtained below can be generalized to any other oligomer 
types). The corresponding model can be expressed as: 

> X -— )• 



X 2 + P f ^ A, 
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where the reactions from the top represent the monomer translation and degradation, dimer- 
ization, dimer degradation, and promoter-binding-unbinding reactions. denotes the dis- 
sociation constant. The time evolution of the total transcription factor (monomer units) 
(Y = X + 2X 2 + 2P b ) is governed by 

dY 

— = a - 7X - 27 2 X 2 . 
at 

The response time constant of Y becomes 



t Pt = - 



dja - 7 X - 2 72 X 2 
dY 



-1 



dX dX 2 

7"^7 + 2 72 



n -1 



dY 



dY 



Under the assumption of the equilibrium of X 2 and P b , we obtain 

M 1^2 1 n PtX% 



Xn 



X\ and P b 



K d + X 9 . 



P T f(X). 



By substituting the first equation in the above to Eq. (14), we obtain 



7 + 



47 2 fciX 
h + 72 



-1 



dY 
dX 



By using Y = X + 2X 2 + 2P b and Eq. fll5[ ), we obtain 

2f{X) 



T~P T 



T 



*!2+72 



p 



T 



where Tq denotes the time constant in the case without any promoter: 



To 



1 , 4fciX 



7 + 



&2+72 



(14) 



(15) 



(16) 



We note that the steady state concentration of X is independent of Pt- Equation (16) shows 
that the extra capacitance C is proportional to Pt (by noting that Tp T = RCp = R{C + C) 
and r = RC). 



B. Degradation-tagged TF 

We consider a monomer TF under directed degradation by proteases. The corresponding 
model can be expressed as: 



X + E XE E 

k 2 [XE] 

K d 

X + P f ^ P b , 
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with the total number of the proteases Et{= E + fixed. The time evolution of the 

total TF (Y = X + [XE] + P b ) is governed by 



dY 
~dT 



a 



l[XE]. 



The response time constant of Y becomes 





~d{a-~f[XE})~ 


-1 


" d[XE]~ 


-l 


Tp T = - 


dY 




7 dY 





Under the assumption of the equilibrium of [XE] and Pb, we can express [XE] and Pb in 



terms of X as like Eq. (15). By following the same way as above, we can prove that the 



time constant increases by the amount proportional to Pt as like Eq. (16). 



C. Multiple operators 

Consider a dimer TF that can bind two operator sites with different affinities. The binding- 
unbinding process between the TF and the promoter can be modeled as 

X + P Q # P x 
X + P ^P 2 

Pi # P12 
P2 ^ P12, 

where Pq, Pi, and Pyi denote promoters that are free, occupied at the z-th operator site, and 
fully occupied, respectively. The equilibrium constant of each reaction from top to bottom is 
denoted by K\, K%, K 3 , and K4. By using that the total promoter concentration Pt equals 
P + Px + P 2 + P m we obtain 

P T X 



P b = P 1 + P 2 + p : 
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X + K d 



P T f(X), 



(17) 



where 



Kd^^ + Ki + K^K-,]- 1 . 
We note that K1/K3 = K2/K4. The response time constant rp T , given by 



Tp 



'dja-jXy 
dY 



T 



dY 
~d~X' 
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is re-written as 

Tp T =T (l + f'(X)P T ). (18) 

by using Y = X + P&. If there were n-operators and the TF is an oligomer, the same 



derivation can be applied resulting in Eq. (18). 
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Figure 1 — Gene circuit modules and their interface process 

(A) The repressilator (Module 1) regulates multiple copies of a downstream module 
(Module 2). The multiple copies can be realized by placing the downstream module in a 
plasmid. TetR repressors (output of Module 1) regulate their specific downstream-module 
promoters (input of Module 2). A module interface process (MIP) defines the collection of 
the processes of ieii?-transcription and translation and its specific downstream-module reg- 
ulation. As the number of the downstream-module promoter (Pt) increases, the amplitude 
and period of the oscillation in the repressilator can be changed (B and C). A repressilator 
model was obtained from the BioModels Database BIOMD0000000012 02]. We modified 
the model to lower the expression levels (by changing translation efficiency to 10 and Km 
to 10 and the maximum transcription rate to 3 per min per cell) and to add promoter 
binding-unbinding reactions for TetR repressors (for the detailed model description, refer to 
the SI). 

Figure 2 — Isolated module output 

Translation-degradation processes for X (A) can be described by a simple reaction process 

(B) with a a translation rate and 7X a degradation rate. These processes can be mapped to 
an RC-circuit with R resistance and C capacitance by V out = X, V in = a/7, and RC = 1/7 

(C) . 

Figure 3 — Module interface process 

The output X of an upstream module regulates the downstream-module X-specific pro- 
moter (A). The translation-degradation processes for X and its promoter regulation can 
be modeled as the reaction process shown in B, where Pf, Pb, and Pt denote the numbers 
of free, bound, and total promoters, respectively. The reaction process is mapped to an 
RC-circuit with an increased capacitance by C, which is shown to be proportional to Pt'- 
C = PtC\ with C\ a proportionality constant. This means that each promoter acts as a 
capacitor with a unit load of capacitance, C\. The total capacitance Ct becomes the sum 
of the capacitance C in the isolated case and the extra capacitance C . 
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Figure 4 — Frequency response of the module interface process shown in Fig. 3A 

An oscillatory signal is applied at V^ n with different frequencies (A). The signal gain is 
defined by the ratio of the oscillation amplitude of the output signal (V out ) to that of the 
input (Vin): g(u) = ^y"^ ' anc ^ can ^ e we ^ approximated by g(u) = + ui 2 / (RCt) 2 
[3T] with Ct the total capacitance. As the number of the promoters that the output (TF) 
signal drives (regulates) increases, the cut-off frequency (u c = 1/RCt) decreases (B). We 
assume that the output signal is desired to be operated within a certain frequency range 
between and 1 hour -1 , i.e., the maximum operating frequency ou max is 1 hour -1 . When 
the cut-off frequency matches the maximum operating frequency, the corresponding number 
of the promoters is defined as the fan-out (C). Parameters of the model: Kd = 1 nM [k on = 
10(l/nM/hour), k off = 10 (1/hour)], 7 = 2(l/hour), a = 20(nM/hour). 



Figure 5 — Module interface processes that the fan-out function Eqs. (10) and (11) 
can be applied to 

(A) An oligomer TF degraded by proteases. (B) A TF can bind two different promoter 
plasmids having different binding affinities and different origins of replication. This can 
be mapped to an RC-circuit with two different capacitances connected in parallel. (C) 
An Oligomer TF can bind multiple operators. (D) Each different TF binds to its specific 
operator without affecting the binding affinity of the other. 



23 



Module Interface 



Module 1 



f r ^ 

J I CI I 


w 


^ s 


1 1 




Irr lad i 












rkll tetR | 

I 1 1 




\ 





Module 2 
(multiple copies) 



B ' 05 

U 
i_ 
CD 
Cl 

V) 
0) 

3 

u 

1 



cd 

E 

3 



C = 

CD 
O 
i_ 
CD 
Q. 

(/) 
OJ 

O 

£ 



CD 
.Q 

E 



100 

75 
50 
25 




100 
75 h 



25 




I 


i i 


Lad 


- P T =0 




TetR 






cl 


\ / 


X I 


1 X t F X 


i 




X X / X, ) 



50 100 150 200 250 300 
t (min) 



1 1 


1 1 1 

Lacl 


- P T =100 


TetR 




cl 


V 7 \ V 
1 1 


i i i 



50 100 150 200 
t (min) 



250 300 



FIG. 1: 



24 




FIG. 3: 

26 



A 

Vi 



/I 

> 


AV ir 




V 




J 






V 



out 



Aav 


r \ 


j 


\J 



B. 



Desired Operating Frequency Range 



1 



10 u - 



3^ 

CD 



10- 1 ■ 



P t =100 
P t =200 
P t =300 
P t =400 



10 



-1 



I 



max 




0)(hour" ) 



Cut-off Frequency 

UJ C for P T = 



lo c for P T = 200 



10 L 



10" 



max 



10 L 



10" 



10" 1 10 u 10 1 

w(hour" 1 ) 

Good 



0) r <CJ 





P T =121 

\ - 





10 L 



10" 



I — 1 — 1 II III 


P T =200 


: 

1 
1 

1 







10"' 10 u 10' 

«(hour" 1 ) 

Fan-out = 121 



10" 1 10 u 10 1 

co(hour" 1 ) 

Bad 



FIG. 4: 



27 




28 



